Effects of large disorder on the Hofstadter butterfly 
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Motivated by the recent experiments on periodically modulated, two dimensional electron systems 
placed in large transversal magnetic fields, we investigate the interplay between the effects of disorder 
and periodic potentials in the integer quantum Hall regime. In particular, we study the case where 
disorder is larger than the periodic modulation, but both are small enough that Landau level mixing 
is negligible. In this limit, the self-consistent Born approximation is inadequate. We carry extensive 
numerical calculations to understand the relevant physics in the lowest Landau level, such as the 
spectrum and nature (localized or extended) of the wave functions. Based on our results, we propose 
a qualitative explanation of the new features uncovered recently in transport measurements. 

PACS numbers: 73.43.Cd 



I. INTRODUCTION 

Two-dimensional electron systems (2DES) placed in a 
uniform perpendicular magnetic field exhibit a rich va- 
riety of phenomena, such as the integer^ and fractional^ 
quantum Hall effects.'^ Another well-studied problem is 
that of a 2DES in a uniform perpendicular magnetic field 
subjected to a periodic potential. Even before the discov- 
ery of the quantum Hall effects, Hofstadter— showed that 
in this case, the electronic bands split into a remarkable 
fractal structure of subbands and gaps, the so-called Hof- 
stadter butterfly. Two "asymptotic" regimes are usually 
considered: (i) if the magnitude of the periodic poten- 
tial is very large compared to the cyclotron energy and 
the Zeeman splitting, then one can use lattice models to 
describe the hopping of electrons between Wannier-like 
states localized at the minima of the periodic potential, 
whereas (ii) if the magnitude of the periodic potential is 
small compared to the cyclotron energy, the periodic po- 
tential lifts the degeneracy of each Landau level. In both 
cases, the resulting butterfly structure is a function only 
of the ratio between the flux = BA of the magnetic 
field through the unit cell of the periodic lattice, and 
the elementary magnetic flux (po — hc/e. Remarkably, if 
(jj/ipo of the first asymptotic case is equal to 4>o/4' of the 
second case, their electronic structures are solutions of 
the same Harper's equation— If the periodic potential is 
comparable to the cyclotron energy. Landau level mixing 
must be taken into account; although Landau levels still 
split into subbands, the structure is no longer universal, 
but depends also on the ratio of the periodic potential 
amplitude and the cyclotron energy— 

Experimentally, the case with a small periodic modu- 
lation can be realized more easily. This is because the 
periodic potential is usually imprinted at some distance 
from the 2DES layer; as a result, its magnitude in the 
2DES is considerably attenuated. The interesting cases 
to study experimentally also correspond to small values 
of 4'/4>{) (of order unity), where the butterfly structure 
shows a small number of subbands separated by large 
gaps, and should therefore be easier to identify. Periodic 



modulations have been created using lithographjiSi^ and 
holographic illuminationiiii The lattice constants of the 
resulting square lattices are of order 100 nm. As a re- 
sult, the condition (/i/^o ~ 2 (for instance) is satisfied for 
B K, 0.8 T. This is a very low value, in the Shubnikov- 
de Haas (SdH) regime, not the high-i? quantum regime. 
Significant Landau level mixing and complications from 
the fact that the Fermi level is inside one of the higher 
Landau levels for such small i?-values make the identifi- 
cation of the Hofstadter structure difhcult. 

Recently, a new method for lateral periodic modulation 
has been developed using a self-organized ordered phase 
of a diblock copolymer deposited on a GaAs/AlGaAs 
heterostructuresli The polymer spheres create a 2D tri- 
angular lattice with a lattice constant of about 39 nm. 
The corresponding unit cell area is almost an order of 
magnitude smaller than those achieved in previous ex- 
periments, implying that the condition 0/</'o ~ 2 is now 
satisfied for very strong magnetic fields, B « 6 T. At such 
high magnetic fields the system is in the strong quantum 
regime, and Landau level mixing can be safely ignored. 
For the experimental 2DES electron concentrations, the 
Fermi level is in the spin-down lowest Landau leveli^i As 
a result, this experimental setup appears more promising 
for the successful observation of the butterfly. 

Nevertheless, one must take into account the disorder 
which is present in the system (without disorder, there is 
no integer Quantum Hall Effect - IQHE - to begin with) . 
If the disorder is very small compared to the periodic po- 
tential amplitude, one expects that the subbands of the 
Hofstadter structure are "smeared" on a scale Ti/t, where 
T is the scattering time, and r — > oo as disorder becomes 
vanishingly small. As a result, the larger gaps in the 
Hofstadter structure should remain open at the positions 
predicted in the absence of disorder, and one expects a 
series of minima in the longitudinal conductivity as the 
Fermi level traverses such gaps. The experiment indeed 
shows a very non-trivial modiflcation of the longitudi- 
nal resistivity, with many peaks and valleys appearing 
in what is (in the absence of the periodic modulation) a 
smooth Lorentz-like peakiii However, the position of the 
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minima in pxx do not track the positions of the main gaps 
in the corresponding Hofstadter butterfly structure. In- 
stead, the data suggests that in this experimental setup, 
disorder is not small, but rather large compared to the es- 
timated amplitude of the periodic potential. This is not a 
consequence of poor samples, since these 2DES have high 
mobilities. It is due to the fact that the periodic modula- 
tion is considerably attenuated in the 2DES, leading to a 
small energy scale for the Hofstadter butterfly spectrum 
as compared to Ti/t. As a result, the Hofstadter structure 
predicted in the absence of disorder is of little use for in- 
terpreting the experimental data. One might expect that 
in this case the periodic potential should have basically 
no effect on the disorder-broadened Landau level. This 
is indeed true for the strongly localized states at the top 
and bottom of the Landau level. However, states in the 
center of the Landau level extend over many unit cells 
of the periodic potential, and, as we demonstrate in the 
following, are non-trivially modified by its presence. 

In this paper, we investigate numerically the behavior 
of a 2DES subject to a perpendicular magnetic field, a 
periodic potential and a disorder potential, under con- 
ditions applicable to the experimental system. The ef- 
fective electron mass in GaAs is 0.067me while the mag- 
netic fields of interest are on the order of 10 T. Under 
these conditions, the cyclotron energy fiujc, of the order 
of 200 K, is the largest energy scale in the problem. The 
Zeeman energy g* for these fields is roughly 3 K, 
but electron interaction effects lead to a considerable en- 
hancement of the spin splitting between the (spin po- 
larized) Landau levels, which has been measured to be 
20 K.^^ The amplitude of the periodic potential's largest 
Fourier components is estimated to be of the order of 
1 K, and the scattering rate from the known zero field 
mobility is estimated to be ?i/r ~ 8 K.~ As a result of 
this ordering of energy scales, we neglect Landau level 
inter-mixing and study non-perturbatively the combined 
effects of a periodic and a large smooth disorder poten- 
tial on the electronic structure of the lowest Landau level. 
Previously, the effects of small disorder on a Hofstadter 
butterfly have been perturbatively investigated using the 
self-consistent Born approximation (SCBA)(i4 and the 
combined effect of white- noise disorder and periodic mod- 
ulation on Hall resistance was studied following the scal- 
ing theory of IQHE.^''^ Our results reveal details of the 
electronic structure not investigated previously. 

The two-lead geometry we consider is schematically 
shown in Fig. ^ the finite 2DES is assumed to have peri- 
odic boundary conditions in the y-direction (along which 
the Hall currents flow) , and is connected to metallic leads 
at the X = —Lx/'i and x — -^-Lx/I edges. In particular, 
in this paper we study the effects of the periodic poten- 
tial on the extended states carrying longitudinal currents 
between the two leads, and identify a number of interest- 
ing properties, in qualitative agreement with simple ar- 
guments provided by a semi-classical picture. Our main 
conclusion is that while the beautiful Hofstadter struc- 
ture is destroyed by large disorder, the system still ex- 
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FIG. 1: The two-lead geometry considered: the finite-size 
2DES has periodic boundary conditions in the y-direction, 
and is attached to metallic leads at the x — ±Lx/2 ends. 



hibits very interesting and non-trivial physics. 

The paper is organized as follows: in Section ^ we 
briefly review the computation of the Hofstadter struc- 
ture for a small-amplitude periodic potential. In Sec- 
tion ^Q] we describe the type of disorder potentials con- 
sidered. Section IIVI describes the numerical methods 
used to analyze the spectrum and the nature of the elec- 
tronic states, with both semi-classical and fully quantum- 
mechanical formalisms. Results are presented in Sec- 
tion while Section IVll contains discussions and a sum- 
mary of our conclusions. 

II. PERIODIC POTENTIAL 

To clarify our notation, we briefly review the problem 
of a free electron of charge — e moving in a 2D plane 
(from now on, the xy-plane, of dimension Lx x Ly) in a 
magnetic field B = Be^ perpendicular to the plane, as 
described by 

1 / e \2 1 
2m V c / 2 

In the Landau gauge A = {0,Bx,0), the eigenfunctions 
of the Schrodinger equation Ti.\n,ky,a) — En,(T\n,ky,a) 
are: 

{r\n,ky,a) = ^=e '""^ X., (1) 

with eigenenergies 

En,a = ?iwc (^n +^ ^ ^gf^sBa. (2) 

Here I — y^hc/eB is the magnetic length, lUc = eB/mc is 
the cyclotron frequency, Hn(x) are the Hermite polyno- 
mials and x+i = (1 0), respectively X-i = (0 1) are the 
eigenspinors of a^: cr^Xo- = crxa- 

The degeneracy of a Landau level is given by the num- 
ber of distinct ky values allowed. Imposing cyclic bound- 
ary conditions in the y-direction, we find 

h = J-j^ (3) 
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where j is an integer. The allowed values for j are found 
from the condition that the electron wave-functions, 
which are centered at positions Xj — Pky — P2TTj/Ly 
[see Eq. (^] are within the boundary along the x-axis, i. 
e. —Lx/2<Xj<Lx/2. It follows that the degeneracy 
of each Landau level is iV = LxLyB /(f>o, with 0o — he/ e. 

Consider now the addition of a periodic potential, with 
a lattice defined by two non-coUinear vectors ai and a2, 
such that V{y) — V {Y+nai+mSi2) for any n,m E Z. The 
periodic potential has non- vanishing Fourier components 
only at the reciprocal lattice vectors g = /igi+fcg2, where 
gi • Sij — 2iT5ij and h, k are integers. Thus: 



(4) 



Further, since V{y) is real, it follows that = Vl^. 

In the absence of Landau level mixing, the Hofstadter 
spectrum for both square** 



Vs{x,y) = 2A 



27r 27r 
cos — X + cos — y 
a a 



(5) 



and triangular- 
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periodic potentials, with nonzero Fourier components 
only for the shortest reciprocal lattice vectors, have been 
studied extensively in the literature4*iSiiLi& The param- 
eter defining the spectrum is the ratio between the flux 
(j) — 'B ■ (ai X a2) of the magnetic field through a unit 
cell and the elementary flux (j)Q. For 4>/(j)o — q/p, where 
p and q are mutually prime integers, the original Landau 
level is split into q sub-bands. 

We would like to emphasize a qualitative difference be- 
tween the two types of potentials: the square potential 
in Eq. Q is particle-hole symmetric, since Vs{x,y) — 
~Vs{x + + As a result, the sign of its amplitude 
is irrelevant. On the other hand, the triangular potential 
does not have this symmetry. With the sign chosen in Eq. 
© and A > Q,Vt has deep local minima at the sites of 
the triangular lattice, whereas the maxima are relatively 
flat and located on a (displaced) honeycomb lattice. As 
a result, the sign of Vt is highly relevant. The second fact 
that must be mentioned is that the choice made in Eqs. 
© and © is rather simple, since it aligns the periodic 
potential with the edges of the sample in a very specific 
way. In general, however, one could consider the case 
where the periodic lattice is rotated by some finite angle 
with respect to the sample edges; study of such cases will 
be discussed in future work. Finally, it may seem that 
this choice of periodic potentials is very restrictive also 
because only the shortest lattice vectors have been kept 
in the Fourier expansion. In fact, the methods we em- 
ploy can be directly used for potentials with more Fourier 
components, but their inclusion leads to no new physics. 



III. DISORDER POTENTIAL 

Real samples always have disorder. The current con- 
sensus is that high-quality GaAs/AlGaAs samples ex- 
hibit a slowly varying, smooth disorder potential. In a 
semi-classical picture, the allowed electron trajectories 
in the presence of such disorder follow its equipoten- 
tial linesiiiS Closed trajectories imply localized electron 
states, while extended trajectories connecting opposite 
edges of the sample are essential for current transport 
through the sample (for more details, see Sec. lIV 

In typical experimental setupsjii dopant Si impurities 
with a concentration of ~ 10"'^^ cm~^ are introduced in 
a thin layer of 6 nm in thickness, located 20 nm above 
the GaAs/AlGaAs interface. Typically, up to 10% of the 
Si atoms are ionized. A small fraction of the ionized 
electrons migrate to the GaAs/AlGaAs interface where 
they form the 2D electron gas. The electrostatic po- 
tential created by the ionized impurities left behind is 
the major source of disorder in the 2DES layer. On the 
length-scale we are interested in, there are 10^ to 10^ such 
ionized Si impurities per jim?' . The resulting disorder po- 
tential must be viewed as a collective effect of the den- 
sity fluctuation of the ionized impurities^*^ rather than a 
simple summation of the Coulomb potential of a few im- 
purities. The electrostatic potential from Si impurities 
is compensated and partially screened by other mobile 
negative charges in the system such as, for example, the 
surface screening effect by mirror charges considered by 
Nixon and Daviesi^ An exact treatment of this problem 
is difficult, since one should consider the spatial corre- 
lation of the ionized impurities iSLSS One model used to 
describe such disorder consists of randomly placed Gaus- 
sian scatterersiS This model captures the main feature 
of a smooth disorder potential and supports classical tra- 
jectories on equipotential contours, but it has no natural 
energy/length scales associated with it. As a result, here 
we choose to also investigate a different model of the dis- 
order, which incorporates the smooth character of the 
Coulomb potential in real space. 

We generate a realization of the disorder potential in 
the following way: positive and negative charges, cor- 
responding to a total concentration of lO'^/im"'^ are 
randomly distributed within a volume [—Lx/2,Lx/2\ x 
[—Ly/2,Ly/2\ X [20nm -|- d, 26nm -I- d] above the elec- 
tron gas which is located in the z — Q plane. Here, 
we choose d = 4 nm as an extra spacer since the elec- 
tronic wave-functions are centered about 3-5 nm below 
the GaAs/AlGaAs interface. Since we are not simu- 
lating single impurities but density fluctuations, these 
charges are not required to be elementary charges. In- 
stead, we use a uniform distribution in the range [— e,e] 
for convenience (a Gaussian distribution would also be 
a valid choice) , and sum up all Coulomb potentials from 
these charges, using the static dielectric constant in GaAs 
e = 12.91.^"' The resulting disorder potential has energy 
and length scales characteristic of the real samples. Typ- 
ical contours for such potentials are shown in Sec. Ivl 
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In an infinite system, in the quantum Hall regime, 
the existence of quantum Hall steps implies the exis- 
tence of critical energies at which the localization length 
diverges.^^ This is the quantum analog of the two dimen- 
sional percolation problem in a smooth random land- 
scape, for which there exists a single critical energy^ 
In the case of potentials with electron-hole symmetry 
{V{r)) — 0, the critical energy lies in the middle of the 
band {Ec = 0), leading to percolating path at half filling. 
For a finite mesoscopic sample, however, not only does 
the percolating path (critical energy) Ef, deviate from 
this value, but in samples without a periodic boundary 
condition one need not have a percolating path traversing 
the system in the desired direction. This arises from the 
fluctuations near the edge of a mesoscopic system with 
free boundary conditions. 

We circumvent such a possibility by adding an extra 
smooth potential V'{x,y) to the impurity-induced dis- 
order potential Vi{x,y), such that the total potential 
V — Vi + V is zero on the opposite edges x — ±Lx/2 of 
the sample where the metallic leads are attached. The 
supplementary contribution V'{x, y) can be thought of as 
simulating the effect of the leads on the disorder poten- 
tial, since the metallic leads hold the potential on each 
edge constant by accumulating extra charges near the 
interface. Therefore, physically we expect that the ex- 
tra potential V' decays exponentially over the screening 
length A inside the sample. This implies: 

, V,{-L,/2,y) + V,{L,l2,y) cosh(a:/A) 
^""'^^ 2 cosh(i,/2A) 

Vd-L,/2,y)-VdL^/2,y) sinh(x/A) 
2 sinh(L:,/2A) 

where A is taken to be 100 nm in our calculation. 

In Fig. 121 we plot the average of Fourier transform 
of the magnitude of the random potential 
versus (Z = |q| for the Coulomb model and the Gaus- 
sian model. The Gaussian model is generated by adding 
100 randomly placed Gaussian scatterers on an area of 
3/Ltm X 3/xm, each contributing Ade~^ , where Ad is 
uniformly distributed in [—2, 2] meV, and d is uniformly 
distributed in [0, 0.2] /^m. F(q) is related to V{y) by 
V{v) = X]q where the summation is over all 

the wavevectors involved in the fast Fourier transforma- 
tion. The Gaussian model has an arbitrary energy scale 
which is fixed by the maximum value of the distribu- 
tion Am. Here Am — 2 meV. As can be seen, V{q) of 
both models are decreasing functions of q. The trends 
of decay are exponential at large q. At small q, the two 
models behave differently. Despite the difference, both 
models lead to the same qualitative results, although, 
as expected, minor quantitative differences are present. 
This shows that the physics we uncover is independent 
of the particular type of slowly- varying disorder potential 
considered, and therefore should be relevant for the real 
samples. 
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FIG. 2: Averaged Fourier amplitudes of two types of disor- 
der potential as a function of wavevector q = |q|. For both 
Coulomb and Gaussian model, V{q)'^ is averaged over 116 dis- 
order realizations. The relation between V{q) and V{r) and 
relevant parameters are discussed in the text. The standard 
deviation, {LxLy)~^ (^J drV'^{r)'^ for the Coulomb model is 
3.2 X lO^^eV^, and 2.1 x 10~^eV^ for the Gaussian model. 



IV. NUMERICAL CALCULATIONS 

In this section we discuss the numerical methods we 
use, including derivations of some relevant formulas. As 
already stated, we focus on the case where the ampli- 
tudes of the periodic and disorder potentials are very 
small compared to the cyclotron energy and the Zeeman 
splitting, and therefore inter-level mixing is ignored. 



A. Semi-classical treatment 

The semi-classical approach is valid— for the integer 
quantum Hall effect in the presence of a slowly varying, 
smooth disorder potential and large magnetic fields (such 
as we consider), so that the magnetic length I which de- 
termines the spatial extent of the electron wave-functions 
is much smaller than the length scale of variation of 
the smooth disorder potential, |Vy(r)| <C huJc/l- Then, 
semi-classically the electron moves along the equipoten- 
tial contours of the disorder potential V{r), in the di- 
rection parallel to WV{r) x B. Since the kinetic energy 
is quenched in the lowest Landau level, the total energy 
of the electron simply equals the value of the disorder 
potential on the equipotential line on which its trajec- 
tory is located. As a result, the density of states in the 
semi-classical approach is directly given by the proba- 
bility distribution for the disorder potential, which can 
be calculated by randomly sampling the potential energy 
and plotting a histogram of the obtained values>i242& 

In Sec. we compare the results obtained within this 
semi-classical approach with fully quantum mechanical 
results. As expected, the agreement is good if only the 



5 



disorder potential is present. However, if the periodic 
modulation is also included, the lattice constant a pro- 
vides a new length-scale which is comparable to the mag- 
netic length I, and the semi-classical picture breaks down. 
Quantum mechanical calculations are absolutely neces- 
sary to quantitatively treat this case. 



B. Quantum Mechanical Treatment 

As shown in SecE] for a finite sample of size Lx x Ly 
at a given magnetic field B, the degeneracy of the unper- 
turbed Landau level is = L^LyB /(po — L^Ly/ {iirP). 
Since the disorder varies very slowly, we need to consider 
systems with L^, Ly ^ I to properly account for its ef- 
fects. As a result, the number of states in a Landau level 
can be as large as 10* in our calculations. Storage of 
the Hamiltonian as a dense matrix requires considerable 
amount of computer memory and its direct diagonaliza- 
tion is prohibitively time-consuming. Sparse matrix di- 
agonalization techniques could be employed, but they are 
less efficient when all eigenvectors are needed, and also 
have some stability issues. 

Here we introduce the numerical methods we use to 
compute densities of states and infer the nature (localized 
or extended) as well as the spatial distribution of the 
wave- functions, while avoiding direct diagonalization. 



1. Matrix elements 

Since inter-level mixing is ignored, the Hilbert sub- 
spaces corresponding to different spin-polarized Landau 
levels do not hybridize. Each Hilbert subspace (n, a) has 
a basis described by Eq. containing A^ orthonormal 
vectors indexed by different ky values. 

In order to compute matrix elements of the total 
Hamiltonian in such a basis, we use the following identity 
derived in Ref. id (notice their different sign convention 
for ky. If 0" 7^ cr', the overlap is zero): 



either case, there are Fourier components corresponding 
to Qy — ±27r/a and qy = 0. Since only basis vectors 
for which the difference ky — ky — qy give non-vanishing 
matrix elements, we must choose the length Ly of the 
sample to be a multiple integer of a, the lattice constant. 

The matrix elements of the disorder potential are com- 
puted in a similar way. We use a grid of dimension A^^; x 
Ny to cover the sample and generate the values of the dis- 
order potential on this grid. Then, fast Fourier transform 
(FFT)Si is used to find the long wavelength components 
of the disorder potential corresponding to the allowed 

I^S- (proper care is 



(7) 



where 



values qx,y = 0,±-^,...,± 

taken to define Fourier components so that Vq = V*^). 
The matrix elements of this discretized disorder poten- 
tial are then computed using Eq. ||7J|. In principle, finer 
grids (increased values for A^a: and A'j,) will improve ac- 
curacy. However, they also result in longer computation 
times, since they add extra matrix elements in the sparse 
matrix, corresponding to large wave-vectors. We have 
verified that a grid size of dimension A^^; = A^j, = 72 is al- 
ready large enough to accurately capture the landscape of 
a 3/im X sample and the computed quantities have 
already converged, with larger grids leading to hardly 
noticeable changes. This procedure is also justified on 
a physical basis. First, the neglected large wave- vector 
components describe very short-range spatial features, 
which are probably not very accurately captured by our 
disorder models to begin with, and which are certainly 
not believed to influence the basic physics. Secondly, this 
procedure insures that the actual disorder potential we 
use is periodic in the ?/-direction, since each Fourier com- 
ponent retained has this property. This is consistent with 
our use of a basis of wave-functions which are periodic 
along y. 

The matrix elements of the Hamiltonian within a 
given Landau level (n, a) are then (n, fcy, (T|7i|n, ky, a) = 
EuM + {n,ky\V\n,k'y), where En,a are given by Eq. (0) 
and the matrix elements of both the periodic and the 
disorder part of the potential V are computed as already 
discussed. This produces a sparse matrix, which is stored 
efficiently in a column compressed format. 



2. Densities of States and Filling Factors 



.ml 

«(q) = ( -j^^ 



qx + iqy 



xe 



9S + ' 

1^(1"' 



\Q), 



m and M the minimum and the 



with g = ip(q2 

maximum of n' and n respectively, and Lm" "'''(Q) the 
associated Laguerre polynomial. When band-mixing is 
neglected n = n' and >Cn,„(q) — e^"^^ Ln{Q)- For the 
first Landau level, Lq(x) — 1. 

Eq. Q gives us the matrix elements for the square 
[Eq. 0] or triangular [Eq. ijHJl] periodic potentials. In 



A quantity that can be computed without direct diag- 
onalization is the filling factor i'n,a-iEp) as a function of 
Fermi energy. The filling factor is defined as: 



E„ 



(8) 



where Q{x) is the Heaviside function and A^ is the total 
number of states in the (ti, cr) Landau level. (Since we 
neglect Landau-level mixing, we can define this quantity 
for individual levels.) The filling factor tells us what frac- 
tion of the states in the given Landau level are occupied 
at r = 0, for a given value of the Fermi energy. This 
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corresponds to the average filling factor measured in ex- 
periment and is also proportional to the integrated total 
(as opposed to local) density of states. 

The filling factor is straightforward to compute if the 
eigenenergies i?ri,a,o- are known. However, we want to 
avoid the time-consuming task of numerical brute force 
diagonalization. The strategy we follow is a generaliza- 
tion to Hermitian matrices of the method used in Ref. 

We restate the problem in the following way: as- 
sume we have a Hermitian matrix of size N x N (no 
Landau level mixing), given by the matrix elements of 
M = Ti. ~ Epl in the basis |n, ky, a) (1 is the unit ma- 
trix). Then, i'n,(j{EF) is proportional to the number of 
negative eigenvalues of the matrix M . We now generate 
the quadratic form M = X]fj=i dCj^ij, and transform 

it into its standard form A4 = X^i^i '^dXiP using the Ja- 
cobian method described below. Here, d^'s are all real 
numbers, and the Xi's are linear combinations of the C^'s. 
This is a similarity transformation which retains the sig- 
nature of the matrix. As a result, even though the num- 
bers di are not eigenvalues of M, the number of negative 
eigenvalues equals the number of negative di values. It 
follows that i'n,c7iEp) is obtained by simply counting the 
number of negative di values for the given Ep. 

The Jacobian method is iterative in nature. First, all 
terms containing and Ci are collected and the needed 
complementary terms are added to form the first total 
square so that Ci and are eliminated from the 

rest of the quadratic form M. The procedure is then 
repeated for all C2 and C2 terms (producing 1^2) etc., un- 
til all N values di are found. Computationally, this can 
be done by scanning the lower or upper triangle of the 
Hermitian matrix M only once. The total number of 
operations is proportional to the number of nonzero ele- 
ments of the matrix, meaning that for a dense matrix it 
scales with N'^ (sparse matrices require much fewer op- 
erations). As a result, this procedure is much faster than 
brute force diagonalization which scales with (for us, 
N ~ 10''). The filling factor i/n,a{E) is a sum of step- 
like functions, with steps located at the eigenvalues. By 
scanning E and identifying the position of these steps we 
can also find the true eigenvalues En_a,a-, with the desired 
accuracy. Finally, the total density of states is given by 
p„a{E)^di^nAE)/dE. 



3. Green's functions: extended vs. localized states 

The advanced/retarded Green's functions are the so- 
lutions of the operator equation 



ihuj-n±i6) G'^'^{uj) = 1 



(9) 



where ^ ^ 0+ . (In practice we use a set of small positive 
numbers, and use the dependence on S to obtain results.) 
If the exact eigenstates and eigenvalues of the total 
Hamiltonian Ti. are known, 7i|n, a, cr) = -En^Q-.o-ln, a, ct). 



(no Landau level mixing), it follows: 
|n, a, cr)(n, a, (t| 



G«'^(co) 



^^G^i(^). (10) 



The exact eigenstates can be expanded in terms of the 



basis states \n,ky,a) as 



,a, cr) ^ y^^Cn,a{ky)\n,ky,cr) 



(11) 



Since the states |n, ky, a) are localized near x = kyP [see 
Eq. (Q)], the coefficients Cn,a{ky) describe the probability 
amplitude for an electron in the state \n,a,cr). Knowl- 
edge of these coefficients allows us to infer whether such 
states are extended or localized in the x-dircction, i.e. 
whether they can carry currents between the leads. 

However, as already stated, we wish to avoid direct 
diagonalization. We can still infer whether the Hamil- 
tonian has extended or localized wave-functions near a 
given energy huj in the following way. We introduce the 
matrix elements: 

G^;^{ky,k';io) = {n,ky,a\G^'^iu;)\n,k'a) 



E 



hu — Er, 



±iS' 



(12) 



If Landau level mixing is neglected, Eq. Q can be rewrit- 
ten in the basis \n,ky,a) as: 



^ ^{huj ± i5)Sk^^k'^ - {n, ky, a\n\n, ky, a) 



xGit:^{k';,k'y;uo) = 5u,,k'^. 



(13) 



We use the popular numerical library SuperLUjSS 
based on LU decomposition and Gaussian reduction algo- 
rithm for sparse matrices, to solve these linear equations. 
Consider now the matrix element G^'^(A:niin, ^max; f^) 
corresponding to the smallest ky = fcmin and the largest 
ky = fcniax values. If all wave-functions with energies 
close to huj are localized in the x-direction, it follows 
that |G^'^(fcniin, ^max; i^)] is a vcry small number, of the 
order e~-^^/^('^\ where ^(w) is the localization length at 
the given energy. On the other hand, we expect to see a 
sharp peak in the value of |G^'^(fcmin, ^max; i^)| if is 
in the vicinity of an extended state eigenvalue, since [see 
Eqs. (|11I12|| ] both Cn,a{k^in) and c„^Q(fcinax) are non- 
vanishing for an extended wave-function with significant 
weight near both the —L^/^ and the Lx/2 edges. More- 
over, the height of this peak scales like 1/5, so by varying 
5 we can easily locate the energies of the extended states. 
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4- Green's functions: local densities of states 

We can also use Green's functions techniques to image 
the local density of states at a given energy E. By def- 
inition (and neglecting Landau level mixing), the local 
density of states in the level (n, cr) is: 



-lra{v\GiAE)\v), 



(14) 



where the second equality follows from Eq. H10() . This 
fmiction traces the contours of probability \<j>n,a,a{j^)\'^ 
for electrons with the given energy E. Its direct compu- 
tation, however, is difficult and very time-consuming. 

For the rest of this subsection, the discussion is re- 
stricted to the Lowest Landau level n = (the value 
of a is irrelevant). We know that in the lowest Landau 
level, electronic wave- functions cannot be localized in any 
direction over a length-scale shorter that the magnetic 
length Z. As a result, it suffices to compute a projected 
local density of states on a grid with I x I (or larger) 
spacings. The projection is made on maximally localized 
wave-function, defined as follows. Let Tq = {xo,yo) be a 
point on the grid. We associate it with a vector: 



ko,yo) = X! \l^y){h\^o>yo), 



(15) 



0,ky,a 



where we use the simplified notation \ ky) = 
for the basis states of the first Landau level (see Eq. (Q) 
and we take 



{ky\xo,yo) 



- + ky{xo+iyo) 



L,, 



(16) 



It is then straightforward to show that 



_ iv-vo} 
772 



2Trl 



(17) 

In other words, |xo,?/o) is an eigenstate of the first Lan- 
dau level strongly peaked at r = rp. (The phase factor is 
due to the proper magnetic translation). Wc then define 
the projected density of states [compare with Eq. H14() ]: 

Pp{xo,yo;E) = -Im(xo,?/o|G''^(^^)|a;o,yo), (18) 

TT 

and use it to study the spatial distribution of the elec- 
tron wave- functions at different energies. Strictly speak- 
ing, the local density of states defined in Eq. H14I) can- 
not be projected exactly on the lowest Landau level, 
because the lowest Landau level does not support a S- 
function ({r\n,ky,a) ^ 0, Vn). However, the coher- 
ent states |a;o,?/o) we select are the maximally spatially- 
localized wave functions in the lowest Landau level, and 



have the added advantage that they can be easily stored 
as sparse vectors, because of their Gaussian profiles [see 
Eq. (|16l) ]. Moreover, in the limit Z (B — + oo) where 
|(r|a;o,yo)| ^ - xo)S{y - j/o), the projected density 
of states pp{xo,yo', E) — > po,o-(r;i?). Therefore, for the 
large B values that we consider here, the projected den- 
sity of states pp should provide a faithful copy of the 
local density of states. 

We compute the projected local density of states fol- 
lowing the method of Ref. Is^. Let Uq be the vector with 
elements {ky\xo,yo) obtained from the representation of 
|a:;o,yo) in the \ky) basis [see Eq. (fT^ ]. and let H be the 
matrix of the Hamiltonian H in the \ky) basis. We gen- 
erate the series of orthonormal vectors Uq, Ui, ... using: 



and for n > 2, 
v„ 

bn-2 



Vl 
Ul 



fl"u„_i, 



Huo, 
Vl - aoUl 



vjvi 



a, 







vJiV„ 



The numbers a„ and 5„ can be shown to be real. We 
do not have a "terminator"^^" to end this recursive series. 
Instead, our procedure ends when the orthonormal set of 
vectors Uo,Ui, ... exhausts a subspace of the lowest Lan- 
dau level containing all states coupled through the dis- 
order and/or periodic potential to the state \xo,yo) (i.e., 
all states that contribute to the projected DOS at this 
point). In the presence of disorder, this usually includes 
the entire lowest Landau level. 

Then, the projected density of states is given by 
Eq. p8|l . where the matrix element of the Green's func- 
tion is the continued fraction: 

{xo,yo\G^{E)\xQ,yo) = 



E — i6 — ao — bl E — i6 — ai — bf [. . .] 

(19) 

Because the Hamiltonian is a sparse matrix, the genera- 
tion of these orthonormal sets and computation of Pp{E) 
for all the grid points is a relatively fast procedure. More- 
over, this computation is ideally suited for parallelization, 
with different grid points assigned to different CPUs. 



V. NUMERICAL RESULTS 

In this section we present numerical results obtained 
using these methods. We have analyzed over 20 different 
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v(E) 

Q g _ (a) Disorder only 



FIG. 3: Profile of tfie disorder potential obtained from our 
Coulomb model on a 3.11/imx2.96/im sample, without (up- 
per panel) and with (lower panel) the V^'(r) correction at the 
X = ±La;/2 edges. The disorder potential varies between 
—3 meV and 3 meV, on a spatial length-scale much larger 
than I — 12.03 nm. The critical region containing extended 
states is in the vicinity of E = 0.06 meV. The contours are 
shown for E =0.0575 meV (dashed), 0.17 meV (thick solid) 
and 0.31 meV (thin solid). These energy values correspond 
to classical filling factors iz—OAl, 0.58 and 0.68 in the upper 
panel and 1^=0.45, 0.56, 0.66 in the lower panel. The differ- 
ence is due to the supplementary smooth potential V' . 



disorder realizations for samples of different sizes, and 
all exhibit the same qualitative physics. Here, we show 
results for several typical samples. The lattice constant 
is always a = 39 nm if periodic potential is present, as 
defined by the experimental systemJi. 

For the first sample, we consider (f>/4>Q — 3/2 {B = 
4.71 T). The magnetic length is I — 12.03 nm, and we 




v(E) 

Qgl (c) A = 0.5meV 



-5 5 
E (raeV) 



FIG. 4: Semi-classical (dashed line) and quantum (solid line) 
filling factors for the disorder potential shown in Fig. but 
different amplitudes of the triangular periodic potential (a) 
A=0- (b) A=0.05meV; (c) A=0.5meV and (d) A^5 meV. As 
expected, agreement exists only in the limit A 0. 



2.964/im. With these choices, the Landau level contains 
N = 10108 states. The disorder potential obtained with 
our scheme described in Sec. Illll is shown in Fig. 13 both 
with and without the correction V'{r). An extended 
equipotential line appears, as expected, ai v ~ 0.5. 

In Figs. 0]and[51we plot the filling factor i'(i?) and the 
corresponding total density of states p{E) as a function of 
E (computation details were given in Sec. lIV B^ . These 
quantities are obtained in the semi-classical limit (dashed 
line) and with the full, quantum-mechanical treatment 
(solid line). Results are shown for 4 different cases: (a) 
only disorder potential and (b, c, d) disorder plus a tri- 
angular periodic potential with amplitudes A = 0.05, 0.5 
and 5 meV, respectively. We only plot a relatively small 
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FIG. 5: Semi-classical (dashed line) and quantum (full line) 
density of states calculated from corresponding filling factors 
in Fig. |1| We show only the center of the disorder-broadened 
lowest Landau level, where the density of states is large. 
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energy interval where the DOS is significant, and ignore 
the asymptotic regions with long tails of localized states. 

While the agreement between the semi-classical and 
quantum-mechanical treatment is excellent in the limit 
A 0, the two methods give more and more different 
results as the periodic potential amplitude is increased. 
As already explained, this is a consequence of the fact 
that the magnetic length I is comparable to the lattice 
constant a, leading to a failure of the semi-classical treat- 
ment when this extra length-scale is introduced. In par- 
ticular, in the case with the largest periodic potential 
[panel (d) of Figs. 0] and |3] we can clearly see the ap- 
pearance of the 3 subbands expected for the Hofstadter 
butterfly at 4>/(f>o = 3/2, although the disorder leads to 
broadened and smooth peaks, and partially fills- in the 
gap between the lower two subbands. This picture [panel 
(d)] is quite similar to the density of states that Ref. Q 
calculated using the self-consistent Born approximation. 
This is expected since the SCBA approach is valid in 
the limit of strong periodic potential with weak disorder. 
However, the SCBA approach is not appropriate in the 
limit of moderate or strong disorder, where the higher 
order terms neglected in SCBA are no longer small. For 
disorder varying on a much longer length-scale than the 
periodic potential, one still expects that locally, on rel- 
ative fiat regions of disorder, the system exhibits the 
Hofstadter-type spectrum. However, these spectra are 
shifted with respect to one another by the different local 
disorder values. If disorder variations are small, then the 
total spectrum shows somewhat shifted subbands with 
partially fiUed-in gaps, but overall the Hofstadter struc- 
ture is still recognizable. On the other hand, for moder- 
ate and large disorder, the detailed structure of the local 
density of states from various flat regions are hidden in 
the total density of states. All one sees are some broad- 
ened, weak peaks and gaps superimposed on a broad, 
continuous density of states. 

We now analyze the nature of the electronic states for 
these configurations. We start with the case which has 
only disorder. In Fig.|Slwe plot |G'''^(A;niin, ^max; E)]"^ as a 
function of the energy E, for different values of S (com- 
putation details were given in Sec. IIVB31) . As already 
discussed, extended states are indicated by large values 
of this quantity, as well as a strong (roughly 1/(5^) de- 
pendence on the value of the small parameter S. 

Figure reveals that as S is reduced, resonant be- 
havior appears in a narrow energy interval E = 0.02 — 
0.36 meV. Panel (a) shows that results corresponding 
to S = IQ-^ eV and 5 = 10"^ eV indeed differ by 
roughly 2 orders of magnitude, with 5 = 10~* eV showing 
sharper resonance peaks. The difference between results 
for i5 = 10~® eV and 5 = 10~^ eV shown in panel (b), is 
no longer so definite. The reason is simply that for such 
small S, the denominator in the Green's function expres- 
sion is usually limited by \E — En^a,iy \ and not by 5 [see 
Eq. H12|) ]. and the dependence on 5 is minimal. Only if 
E is such that \E — -E„,a,<T| < ^ can we expect to see a 
5 dependence, and indeed this is observed at some ener- 



gies. Finally, in panel (c) we show the comparison with a 
larger energy interval. The value of the Green's function 
decreases exponentially fast on both sides of the critical 
region, indicating strongly localized states. Here, data 
for 5 — 10^^ eV is a smooth curve, whose magnitude is 
much less than that of the other three values even for lo- 
calized states. This is due to the fact that this 5 is larger 
than typical level spacings. As a result, several levels con- 
tribute significantly to Green's function at each E value, 
and the destructive interference of the random phases 
of different eigenfunctions lead to the supplementary 5- 
dependence. We conclude that the disorder potential has 
a critical energy regime of approximately 0.3 meV width, 
covering less than 5% (in energy) and 20% (in number of 
states) of the disorder-broadened band with total width 
~ 6 meV. The position of the critical energy interval is in 
agreement with the semi-classical results which suggest 
an extended state in the vicinity of E — 0.06meV. By 
comparison with Fig. ^ we can also see that this criti- 
cal regime corresponds to a roughly half-filled band, in 
agreement with the experiment. 

The effect of an additional triangular periodic poten- 
tial is shown in Fig. [7| where we plot the same quantity 
shown in Fig. |H| for a fixed 6 = 10~^ eV and different 
amplitudes ^ = 0, 0.05, 0.5 and 5 meV, respectively. 
These results correspond to a different Coulomb disorder 
potential (not shown), as can be seen from the different 
location of its extended states. Here we see how the nar- 
row critical interval of extended states grows gradually 
as the amplitude of periodic potential is increased and 
finally exhibits the three well-separated extended sub- 
bands expected for (f)/(j)o = 3/2 in the limit of vanishing 
disorder. The three subbands can already be resolved for 
the moderate case A — 0.5 meV, although they are very 
wide and exhibit significant overlap. 

Qualitatively similar behavior is obtained if we use 
the Gaussian scatterers model for disorder. A typical 
realization of this disorder is shown in Fig. |S1 Results 
for the Green's function's values with such disorder are 
shown in Fig. El for cases with pure disorder, and also 
cases with either a triangular or a square periodic po- 
tential. The magnetic field has been doubled, such that 
4'/4'o = 3. Similar to the case shown in Fig. [7| the peri- 
odic potential leads to a widening of the critical regime. 
For large periodic potentials, the expected Hofstadter- 
like three-subband structure emerges again. We conclude 
that Coulomb and Gaussian disorder models show qual- 
itatively similar behavior. 

We now analyze the projected local density of states 
Pp{E) discussed in Sec. IIV B 41 in order to understand 
the reason for this substantial widening of the critical 
region by even small periodic potentials. We consider a 
smaller sample, of size approximately 1.6/imx 1.6/im, and 
compute the projected density of states for 500 equally- 
spaced energy values, on a 60 x 60 square grid and for a 
value S = 10~^ eV. This i5 value is comparable or smaller 
than the level spacing, so we expect to see sharp res- 
onances from the contribution of individual eigenfunc- 
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FIG. 6: Semi-log plot of the amplitude of Green's function 
matrix element between the two edge states near x = ±Li/2, 
as a function of energy. Only the disorder potential of 
Fig. El is present, (a) comparison between S = 10~^ and 
5 = 10~* results; (b) comparison between 5 = 10~* and 
5 = 10~^ results; (c) comparison between results correspond- 
ing to 5 = 5 = 10"^10"^10"® and 10"®. (the last three 
curves are indistinguishable to the eye on this scale.) All 5 
values are in eV units. 
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FIG. 7: The effect of a triangular periodic potential on the 
critical energy regime. The disorder potential used here (not 
shown) supports a narrow interval of extended states centered 
at about —0.6 meV. As the amplitude A of the periodic po- 
tential increases, the range of extended states increases dra- 
matically. The left panel shows results for disorder-only and 
two relatively weak periodic potentials, while the right panel 
shows two larger periodic potentials where the three-subband 
structure expected for 4'/(j>o ~ 3/2 is clearly seen. 



tions as we scan the energy spectrum. Each computation 
generates a large amount of data (roughly 24M), corre- 
sponding to the 500 plots of the local density of states 
at the 500 values of E. Since we cannot show all this 
data, we select a couple of representative cases and some 
statistical data to interpret the overall results. 




FIG. 8: A disorder potential of Gaussian type on a roughly 
3/im X 3pm square. The three lines are equipotential con- 
tours close to the critical regime, with energies of -0.1 meV 
(dashed), meV (thick solid) and 0.1 meV (thin solid). Cyclic 
boundary condition are applied in the y direction. 
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FIG. 9: Green's functions for a sample with Gaussian disor- 
der and various periodic potentials. The calculation included 
20216 states with (j)/4>o = 3. Similar to results shown in Fig.|5] 
we see that the periodic potentials widen the critical region. 



Figures^land^Jshow some of our typical results. The 
two figures are calculated for the same Coulomb-disorder 
potential, for values of i? = —0.504 meV (at the bottom 
of the band) and E — —0.124 meV (close to, but be- 
low the band center) respectively. Each figure contains 
four panels, panel (a) shows the profile of the disorder 
potential as well as an equipotential line (solid black) 
corresponding to the value E considered; the other three 
panels show the projected density of states pp{E) for (b) 
pure disorder; (c) disorder plus triangular periodic poten- 
tial with A ~ Q.l meV; (d) disorder plus square periodic 
potential with A = 0.1 meV. In Fig.^| this equipoten- 
tial line (which traces the semi-classical trajectory of elec- 
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FIG. 10: (color online) Projected local density of states 
pp{r;E) for E — —0.504 meV. Panel (a) shows the profile of 
the disorder potential, and the equipotential contour (black 
line) corresponding to iJ = —0.504 eV. The other three panels 
show pp{r;E) for (b) disorder only; (c) disorder plus trian- 
gular periodic potential with ^4 = 0.1 meV; (d) disorder plus 
square periodic potential with A = 0.1 meV. The width and 
length of the sample are both 1.6 fim, and (?!)/</>o ~ 3/2. In- 
creased brightness corresponds to larger values. 



trons with the same energy E) surrounds local minima 
of the disorder potential, suggesting localized electron 
states at such low energies. Indeed, this is what panels 
(b), (c) and (d) show. The projected density of states 
Pp{E) is large (bright color) at the positions where elec- 
trons of energy E are found with large probabilities. For 
pure disorder, we observe only closed trajectories (local- 
ized states), whose shape is in excellent agreement with 
the semi-classical trajectory, as expected. If a moderate 
periodic potential is added, the wave- functions spread 
over a larger area, and nearby contours sometimes merge 
together. Instead of sharp lines, as seen in panel (b), the 
contours now show clear evidence of interference effects of 
the wave- functions on the periodic potential decorating 
the electron reservoirs. Some periodic modulations can 
also be observed in the background of panels (c) and (d), 
especially for the square potential. These are not the di- 
rect oscillations of the periodic potentials, since the grid 
we use to compute these figures has a linear size equal to 
7/10 of the period a = 39 nm of the periodic potential. 
Capturing detailed behavior inside each unit cell would 
require a much smaller grid, which is not only time con- 
suming, but also violates the requirement that the grid 
size be of order I or larger. 

Figure lTTl for an energy close to the band center, shows 
the same characteristics. For pure disorder, the elec- 
trons at this energy trace a sharp contour very similar 




FIG. 11: (color online) The same as in Fig. 1101 but for an 
energy E = —0.124 meV close to the band center. 

to the corresponding equipotential line shown in panel 
(a) . Electrons are still not delocalized, since this contour 
does not connect either pair of opposite edges. However, 
addition of the periodic potential now leads to extended 
states for both types of periodic potentials [(c) and (d)] 
at roughly E — —0.124 meV, demonstrating the widen- 
ing of the critical region with the addition of a periodic 
potential. 

Physically, one can understand this spread of the wave 
function in the presence of the periodic potential using 
the semi-classical picture>ii If only a smooth disorder 
potential is present, the equipotential at any energy E 
must be a smooth, continuous line. However, if a pe- 
riodic potential with minima —Vm and maxima Vm is 
superimposed over disorder, the new equipotential line 
now breaks into a series of small "bubbles" surrounding 
the disorder-only contour. This happens throughout the 
area defined by the equipotentials E — Vm and E + Vm 
of the disorder potential, since the addition of the peri- 
odic potential leads new regions in this area to have a 
total energy E. Quantum mechanically, we expect some 
tunneling inside this wider area and this is indeed what 
we observe in Figs. 1101 and II II This mechanism suggests 
enhanced delocalization on both sides of the critical re- 
gion as localized wave functions spread out over larger 
areas, as well as a widening of the critical region itself, 
in agreement with our numerical results. 

This spreading of the wave functions in the presence 
of the periodic potential can also be characterized by 
counting, at a given energy E, the number of grid points 
r which have a value pp(r;E) > pc, where pc is some 
threshold value. For sufficiently large pc, this procedure 
counts grid points where electrons with energy E are 
found with large probabilities, thus, in effect it character- 
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FIG. 12: Number of grid points for which pp(r, E) > 100 as 
a function of the energy E. This quantity characterizes the 
"spatial extent" of the wave function. The upper panel com- 
pares results for disorder only (A=0) and disorder plus trian- 
gular potentials with A ± 0.01 meV. The difference observed 
for the two signs is a consequence of particle-hole asymmetry 
of the triangular potential. The lower panel shows results for 
disorder only and disorder plus a square potential. 



izes the "spatial extent" of the wave functions. The re- 
sults of such counting are shown in Fig.^lfor 500 energy 
values corresponding to the disorder potential analyzed 
in Figs. Uniandini There are a total of 60 x 60 = 3600 
grid points on the sample. For the case of pure disorder 
(black line) we see that the largest values are found at en- 
ergies just below 0, where the extended states (the critical 
region) are found for this particular realization of disor- 
der. Because it is a smooth, sharp line, even the most 
extended trajectory has significant probabilities at only 
about 10% of the grid sites. For both higher and lower 
energies, this number decreases very fast, indicating wave 
functions localized more and more about maxima or min- 
ima of the disorder potential, as expected. Addition of a 
small periodic potential increases this number substan- 
tially, clearly showing the supplementary spreading of the 



wave functions in the presence of the periodic potential. 

Figure E| shows this effect for three types of periodic 
potential: triangular lattices with A > and A < (up- 
per panel), and square lattice in the lower panel. All 
three cases show significant enhancement, as compared 
to the pure disorder case. In addition, we see that while 
the square potential gives a fairly symmetrical enhance- 
ment, the triangular potential does not, with curves for 
zLA not overlapping. This is a consequence of the asym- 
metric shape of the periodic potential, which has differ- 
ent values for its minima and maxima \Vm\ 7^ |^Af|j 
well as different arrangements for the points where min- 
ima/maxima appear (triangular lattice vs. honeycomb 
lattice). Fig. E| clearly shows that A > favors in- 
creased delocalization below the critical energy regime, 
while A < favors increased delocalization above it. 

The reason for this different response to the two signs 
of the triangular potential can be nicely explained within 
the semi-classical framework. In Fig. E| we show the 
equipotential lines corresponding to filling factors v = 0.3 
(well below critical region) and v — 0.7 (well above the 
critical region) for a realization of Coulomb disorder (not 
shown) plus a triangular potential with A > 0. Areas 
with energy below the equipotential value are shaded. In 
this case we can clearly see that instead of the continu- 
ous, smooth trajectory expected for disorder-only cases, 
there are also extra "bubbles" regions connecting the ar- 
eas between such contours. Since the choice A > leads 
to deep minima at —Vm = —6A with triangular arrange- 
ment and relatively flat maxima at +Vm = 3^4 with hon- 
eycomb arrangement [see Eq. ©], it follows that the tri- 
angular (honeycomb) "bubbles" region appear roughly in 
the area bounded by the equipotentials E and E+Vm (re- 
spectively, E — Vm and E) of the pure disorder potential. 
At low filling factors, the pure disorder E equipotential is 
a collection of closed contours surrounding local minima 
[see panel (d) of Fig. ^| for an illustration] . It follows 
that for the choice Vm > Vm, the more extended region 
with triangular "bubbles" will be found outside these "is- 
lands" and will lead to a spread of the wave function over 
considerably larger areas, as indeed seen in the upper 
panel of Fig. E| On the other hand, at large filling fac- 
tors the disorder-only contours are "islands" surrounding 
the maxima of the disorder potential. In this case, con- 
tours between E and E + Vm are inside the E contour, so 
the triangular "bubbles" region does not help to connect 
various "islands" as before. The honeycomb "bubbles" 
regions does this, but because Vm < Vm the extension 
of the wave function between "islands" is significantly 
smaller in this case. 

In the quantum-mechanical case one expects interfer- 
ence (due to tunneling) effects between the small "bub- 
bles" regions, and therefore a wave function which is ex- 
tended over their entire indeed we observe to 
be the case in Figs. 1101 and 1111 In other words, one ex- 
pects that a triangular potential with A > will lead to 
considerable increase of the localization length, and re- 
spectively widening of the critical energy region, at filling 
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FIG. 13: Equipotential contours at filling factor 1/ — 0.3 
(upper panel) and v = 0.7 (lower panel) for a 3fim x 3fj,m 
sample with disorder plus a small triangular periodic poten- 
tial with ^ > 0. The shaded regions correspond to energies 
below the respective equipotential. In the semi-classical ap- 
proximation, the shaded areas are filled with electrons, with 
the maximum density of l/27r/^, whereas the white areas are 
completely devoid of electrons. Quantum- mechanically, one 
expects interference in the regions with small periodic "bub- 
bles", induced by the periodic potential (see Figs. 1101 and 
[IB- 



factors below one-half, whereas A < Q wiU favor delocal- 
ization at filhng factors above one-half, as seen in Fig. 1121 
This asymmetry is therefore clearly a consequence of the 
asymmetry of the triangular potential, and is absent for 
a square potential with only lowest order Fourier coef- 
ficients, which possesses electron-hole symmetry. This 
should have clear implications for the transport proper- 
ties of the system. 



VI. SUMMARY AND CONCLUSIONS 

In this study we have investigated the effects of 
moderate-to-large smooth disorder on the Hofstadter 
butterfly expected for 2DES in a perpendicular magnetic 
field and a pure periodic modulation. The parameters of 
our study are chosen so as to be suitable for the interpre- 
tation of recent experiments on a two dimensional elec- 
tron system in GaAs/AlGaAs heterostructure with peri- 
odic modulation provided by a diblock copolymer, The 
experiment shows that (i) the longitudinal resistance Rxx 
is still peaked approximately at half filling; (ii) there are 
many reproducible oscillations in Rxx, indicating non- 
trivial electronic structures in the patterned sample; (iii) 
the distribution of these oscillatory features is asymmet- 
ric, with most of them appearing on the high magnetic 
fields (i.e. low filling factors ly < 0.5) side of the peak 
of Rxx ; and (iv) the temperature dependence of Rxx in- 
dicates that the asymmetric off-peak resistance is ther- 
mally excited, whereas the central Rxx peak (close to half 
filling) has metallic behavior. 

These observations cannot be explained on the basis of 
the Hofstadter structureiii This is not surprising, since 
one expects that large disorder will modify this struc- 
ture considerably. Effects of small disorder on the Hof- 
stadter butterfly had been investigated previously using 
SCBA,^^ but this basically perturbational approach is 
not appropriate for the case of moderate-to-large smooth 
disorder. Instead, we identify and use a number of tech- 
niques which give the exact solution (if electron-electron 
interactions, as well as inelastic scattering are neglected) 
while avoiding brute force numerical diagonalizations. 

Our results demonstrate that while the Hofstadter but- 
terfly is destroyed by large disorder, the effects of the 
periodic potential are non-trivial for states near the crit- 
ical regime. Firstly, they lead to a significant increase in 
localization lengths of the localized states at mesoscopic 
(/im) length scale and induce an effective widening of the 
critical regime near the critical regime. This is achieved 
through a spreading of the electron wave-function on 
the flat regions of the slowly varying disorder potential, 
where their behavior is dominated by the periodic mod- 
ulation. This regime shows an interesting transition be- 
tween the pure disorder and the pure periodic potential 
cases. In the case of pure disorder, the semi-classical 
approach tells us that at finite filling factors, some ar- 
eas of the sample are fully occupied by electrons with 
the maximum possible density of l/(27r^^) (these are the 
areas where the disorder potential has minima) whereas 
other areas are fully devoid of electrons (areas where the 
potential has maxima) and the boundary between such 
regions is very sharp. On the other hand, for a pure pe- 
riodic modulation all wave functions have translational 
invariance with the proper symmetry, and therefore elec- 
tron densities are uniform over the entire sample (up to 
small periodic modulations inside each unit cell). 

When both types of potential are present, with disor- 
der being dominant, our results show three types of areas. 
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There are regions which are fully occupied and regions 
which are completely devoid of electrons, as in the case 
of pure disorder. However, the periodic potential leads 
to a widening of the boundary between the two, where 
the wave functions interact with several oscillations of 
the periodic modulation and therefore have some partial 
local filling. As the critical regime containing wave func- 
tions percolating throughout the sample is approached, 
this spreading of the wave function becomes dominant in 
establishing the transport properties of the system. 

An equivalent way to say this is that the main effect of 
the periodic potential is to provide bridging between the 
fully-occupied electron "puddles" created by the disorder 
potential. Since the connecting areas are relatively flat, 
the local wave functions respond to the local periodic po- 
tential, and therefore locally have a Hofstadter-butterfly 
like structure. If the partial filling factor in such a region 
is inside the gap of the local Hofstadter butterfly struc- 
ture, one expects no transport through this local area. 
This should result in a dip in the longitudinal transport, 
since in such cases the periodic potential will not trans- 
port electrons from one "puddle" to another one. By 
contrast, if the local filling factor in such a region is in- 
side a subband of a local Hofstadter structure, this area 
will establish a link between different "puddles" and thus 
help enhance the transport through the sample. Trans- 
port in this regime should show strong thermal activated 
behavior, in contrast to metallic transport in the critical 
regime where the wave functions connect opposite edges 
of the sample. 

As a result, one expects a series of local minima and 
maxima in the longitudinal resistivity on either side of 
the central peak induced by the extended states (criti- 
cal regime). Furthermore, for an asymmetric triangular 
potential, this response should be strongly asymmetric, 
with the effect most visible on one side of the central 
peak. (One must keep in mind that since tunneling leads 
to exponential dependencies, even small differences in the 
extent of the wave functions can have rather large effects 
on pxx)- Such an asymmetry should also be present in 
longitudinal conductance of finite but low temperature, 
e. g. in the hopping regime which is sensitively depen- 
dent on the nature of the localized wavefunctions, as is 
indeed seen experimentallyjii 

To summarize, our qualitative explanation for the var- 
ious experimental features are as follows: 

(i) The Rxx peak is roughly at the center of the band 
because the weak periodic potential cannot establish a 
Hofstadter-like structure over the whole band. Instead, 
low and high ly states are strongly localized and do not 
transport longitudinal currents. 

(ii) New extended states induced by the periodic poten- 
tial are responsible for the reproducible peaks and valleys 
appearing in Rxx- 

(iii) The periodic potential also leads to the expansion 
of localized wave functions, which contribute to the ther- 



mally activated conduction at lower filling factors. The 
detailed structure of the wave functions gives rise to the 
oscillations of the off-peak Rxx, similar to conductance 
fluctuations.'^^ Finally, 

(iv) the asymmetry in Rxx is a manifestation of the 
asymmetry of the triangular potential, which has a 
stronger effect at low filling factors than at high filling 
factors for A > 0. We predict that this asymmetry should 
be absent for a symmetric square periodic potential. 

The weak point in our calculation is that we are un- 
able to accurately model the potential in the real samples, 
because various screening effects have not been properly 
taken into account. Also, we have no quantitative in- 
formation about the magnitude of the periodic poten- 
tial in the 2DES layer, because of the additional strain^ 
contribution induced by the periodic decoration. As a 
result, we only claim qualitative agreement with the ex- 
periment, although our investigations show the same type 
of behavior for various types of disorder potentials and 
various (small-to-moderate) strengths of the periodic po- 
tential. The most direct check of this work would be 
an experimental demonstration that thermally activated 
conduction appears symmetrically on both sides of the 
Rxx peak for a periodic potential with square symmetry 
and primarily lowest Fourier coefficients. 

Limited computer resources restrict our calculations to 
samples no larger than 3/zmx3/xm, while the sample used 
in experiment has a size of 20/zmx20/zm. From a theo- 
retical point of view, it is interesting to ask what is the 
thermodynamic limit. For pure disorder, it is believed 
that in this limit the typical size of wavefunction diverges 
at a single critical energy. Since we cannot pursue size- 
dependent analysis for samples larger than 3/imx3/im, 
we do not know whether the small periodic potential will 
lead to a finite size critical regime in the thermodynamic 
limit, although this seems likely. From an experimen- 
tal point of view, the interesting question is whether the 
Hofstadter structure can be observed at all. Our studies 
suggest that this may be possible for small mesoscopic 
samples, where the slowly- varying disorder has less ef- 
fect. Alternatively, one must find a way to boost the 
strength of the periodic modulations inside the 2DES. 
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